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Abstract 

A first principles approach is presented for calculations of optical - ultraviolet (UV) spectra 
including excitonic effects. The approach is based on Bethe-Salpeter equation calculations using the 
NBSE code combined with ground-state density- functional theory calculations from the electronic 
structure code ABINIT. Test calculations for bulk Si are presented, and the approach is illustrated 
with calculations of the optical spectra and birefringence of a-phase SiC>2 and the rutile and anatase 
phases of TiC>2- An interpretation of the strong birefringence in TiC>2 is presented. 

Ch ' PACS numbers: PACS 

O' 

o 

> 

o 

(N 

o> 
o 

00 
O 

> 

X 

S3 



I. INTRODUCTION 

As in many subfields in condensed matter physics, the application of density-functional 
theory (DFT)^^ has significantly influenced studies of material optical properties. How- 
ever, treatments based on Kohn-Sham dipole transitions alone only provide a beginning for 
modeling optical spectra, since a ground-state one-electron approach cannot account for a 
range of excited-state and optical effects. Thus quantitative calculations of optical spectra 
require the integration of DFT with a number of many-body techniques. Local field cor- 
rections from the random phase approximation (RPA) alone are inadequate to reproduce 
near-gap optical spectra, because they do not take into account self-energy and excitonic ef- 
fects. Consequently a quantitative theory and interpretation must also include quasi-particle 
self-energy effects based e.g., on the GW- method,- 1 ^ and particle- hole interactions based on 
the Bethe-Salpeter equation (BSE).- The success of the BSE method has effectively settled 
the long-running discussion of the origin of many low energy spectral features.-^ For ex- 
ample, for the rutile phase of TiC>2, early studies could only speculate about the relative 
importance of excitonic effects versus band structure topology in the dominant threshold 
absorption,- whereas we show here that these features are dominated by excitonic effects. 

While the BSE is often necessary for accurate calculations of optical spectra, the ap- 
proach has had limited applications, largely due to a lack of efficient computational meth- 
ods for general systems. To address this need, we have developed an approach based on 
an interface, here dubbed AI2NBSE, between the BSE code NBSE developed at National 
Institute of Standards and Technology (NIST) 10 and the general purpose DFT electronic 
structure code ABINIT.— One of the objectives of this work is to improve the availability 
of BSE codes. Our interface also provides a comparison to the BSE codes EXG^, and 
EXC1ITING, 13 " and other recent codes,— and also to approaches which model electron-hole 
interactions within time-dependent density- functional theory (TDDFT). 15 The latter cal- 
culations are usually simpler than BSE ones, but in practice are limited by an incomplete 
knowledge of the exchange-correlation functional and the neglect of damping effects. The 
overall strategy of our interface differs from these approaches in several respects. For ex- 
ample our AI2NBSE interface achieves efficiency in the BSE calculations through the use 
of the Hybertsen-Levine-Louie^ 1 ^ dielectric screening, and can also treat finite momentum 
transfer. Also, the interface requires only generic input, and thus can be adapted to other 



ground-state and BSE codes. 

The interface is tested on bulk Si, yielding results in good agreement with other ap- 
proaches. In particular we find that the calculations based on ABINIT are in excellent 
agreement with those from the optimal basis function (OBF) code.— As initial applica- 
tions, we report optical spectra and anisotropic optical properties of the common rutile 
and anatase phases of Ti02, and for the a-phase of SiC>2. These are important materials 
for many applications.— 1 ^ However, we are not aware of earlier calculations for TiC>2 that 
include electron-hole interactions, which are needed to reproduce their rich dependence on 
polarization and phase. 

In the remainder of this paper we briefly summarize the key formulas describing excitonic 
effects and other features of NBSE. We then briefly describe our integration of ABINIT and 
NBSE, using bulk Si as an illustrative test case, with some further details in an Appendix. 
Finally, calculations are presented for the birefringent spectra of the a-phase Si0 2 , and 
similarly for the rutile and anatase phases of TiC>2, followed by a summary and conclusions. 

II. BSE AND EXCITONIC EFFECTS 

A. The Exciton Secular Equation 

In this paper we only briefly summarize the BSE formalism, following the notation and 
theoretical developments of Shirley et al— Formally the BSE provides a complete theory 
for optical spectra through a hierarchy of equations derived from the two-particle Green's 
function.— 122 However with certain approximations to the electron-hole interaction, the BSE 
can be reduced to an eigenvalue problem of an effective particle-hole Hamiltonian^ 

H\f} = [H le + H 2e }\f} = n f \f), (1) 

where the eigenstates |/) are given by a superposition of particle-hole basis states^ 
|rm'k(q)), i.e., 

I/) = ^2*l>f(nn'k(q)) |nra'k(q)) 

nn'k 

|nn'k(q)} = a nk oi, k+? |0), (2) 

and throughout this paper we use atomic units (e=h=m=l) . Here if)f(nn''k(q)) is the 
amplitude of a given particle-hole (or excitonic) state with Bloch crystal momentum index 



k and momentum transfer q, the index n runs over all occupied valence bands, the index 
n! runs over unoccupied bands, and |0) denotes the many-body ground-state with energy 
Qq = 0. For optical absorption q is usually negligible. However for inelastic x-ray scattering 
and other spectroscopies and for computational purposes, it is desirable to retain the explicit 
momentum-transfer dependence. The single-particle contribution to the Hamiltonian for a 
particle-hole pair is diagonal in the quasi-particle/hole basis, so that one has 

H le \nn'k(q)) = (£ n 'k+q - Kk) |wi'k(q)), (3) 

where the quasi-particle energies E nk are Kohn-Sham eigenvalues e n k plus quasiparticle 
self-energy corrections 

E n k = e„k + S n k. (4) 

Here £ ra k is the self-energy calculated in the GW approximation for which efficient approx- 
imations have been developed.— 1 ^ 

The electron-hole interaction contribution includes both "direct" and "exchange" cou- 
plings H 2e = V D + V X , 

H 2e \nn'k{q)) = J^ [V D {nn'k,n"n'"k!; q) 

n"n'"k 

+ V x (nn'k, n"n'"k'- q)] |nV"k'(q)>, (5) 

with matrix elements defined explicitly below. Once the BSE secular equation is solved, 
the optical constants may be obtained formally using a Fermi golden rule expression in 
terms of the excitonic final states coupled to the current-like operator J' . In the NBSE 
code, however, these properties are calculated using resolvent techniques. In particular the 
imaginary part of the dielectric tensor is given in terms of resolvents 

Ime^H = -4vrlm [(0| J> - H + ir}]' 1 J' v \0) 

-{OlJU-u-H-ir,}- 1 ^)]. (6) 

In terms of the particle-hole states, the current-like operator coupling to the ground-state 
is given approximately by 

j;io)-^|Wk(q))^±^^, (7) 

nn'k ^ 



where J M is the /x-th component of the current operator and ip n \<_ are approximated as Kohn- 
Sham one-particle states. For small q, the matrix elements of J^ are approximated by 

-(VVk+q|<//J^nk) 
10 

£w/k+q ~ £nk ) (^k +q |e iq - r |Vn k ) 

~-(VVk +q |e iq - r |^k). (8) 

where q = qfi. In the NBSE code these resolvents in the above expressions for the dielectric 
tensor are calculated using an efficient iterative Lanczos algorithm.— 



B. Electron-hole Interaction 

The interaction kernel H 2e of Eq. (5) accounts for two processes which scatter an electron 
from band n' to band n" f , and a hole from band n to band n" . The first is the attractive 
direct screened Coulomb interaction between the electron and the hole, and the second is 
the repulsive unscreened exchange interaction. Each of these contributions can be written 
in terms of two-particle integrals between electron and hole Kohn-Sham orbitals ^> n k(x) and 
■0n'k+q( x )- Matrix elements of the direct term are given by 

V D {nn'k,n"n'"k';q) « - / d 3 x^*,,, k , +q (x)Vv k + q ( x ) 

x |rfVC k (x')t„ k ,(x')f(x,x'; U = 0). (9) 

Here the electron-hole interaction W{k,x';u) is the screened Coulomb attraction as men- 
tioned above, i.e., W = e _1 (dj)l / , which in an exact theory should include the frequency 
dependence of the dielectric response. However the NBSE code approximates the screening 
with the static, spatially dependent Hybertsen-Levine-Louie dielectric function.— 1 ^ This 
model maps the local density in real solids to the density dependence in a homogeneous 
system. These calculations require the ground-state density, whose Fourier coefficients are 

-°( G ) = Jf E /nkV>„k(G + G')Ck(G'), (10) 

nk,G' 

where J„k = 6(sf — £nk) is the occupation of state rik, N is the number of unit cells in 
the ABINIT calculation, which is the same as the number of k points in the calculation, 



and spin degeneracy is assumed. Here ^nk(G) are the Fourier coefficients of the Bloch wave 
expansion 17 in reciprocal lattice vectors 

^ k (x) = e ik ' x J2 ^k(G) e iGx . (11) 

G 

Because the dielectric function is modeled locally, the exciton amplitudes must be trans- 
formed from the Bloch basis to a local basis with coordinates x and x'. These local coor- 
dinates can be considered as transform analogs of the band and wave-vector indices n, n', 
and k.— 

Similarly matrix elements of the exchange term are 

Vx{nn'k, n"n'"k';q) « +2 / d 3 x -^*,,, k , +q (x)?/y, k ,(x) 

xj dV^xOVnWxOllx-x'l. (12) 

As noted, for example by Hybertsen and Louie, the exchange term in the BSE should be 
unscreened. 

The matrix dimension of the electron-hole interaction Hi e is generally very large. For 
example, for the calculations for SiC>2 presented below, 216 k points are sampled, and there 
are 24 doubly degenerate valence (hole) states and 26 conduction (electron) states. For this 
case, there are over 10 5 electron-hole pairs, each requiring a representation in both periodic 
and local bases. 

C. AI2NBSE Interface 

The above theory has been implemented in a modular code which uses the output 
of ground-state density-functional theory calculations from the electronic structure code 
ABINIT as input to Bethe-Salpeter equation calculations using the NBSE code. The in- 
terface serves as a driver for both ABINIT and NBSE, starting from a single input file, 
and also constructs the various physical quantities and arrays needed in the calculations. No 
explicit changes in the structure or coding of either ABINIT or NBSE are used. Both the 
interface and documentation are available from the authors.— Additional details are given 
in the Appendix. 
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FIG. 1: Imaginary part of the dielectric function for Si using NBSE and two different ground-state 
inputs; one set of inputs is from OBF (dashed line) and the other from ABINIT (solid line). For 
comparison the experimental spectrum is also plotted (crosses). 

III. OPTICAL SPECTRA OF SILICON 

As a quantitative test, illustrative results from AI2NBSE for bulk Si are presented in Fig. 
1. The two spectra compared are each calculations of the imaginary part of the dielectric 
function for bulk Si using NBSE: in one case the ground-state quantities above are calculated 
with the optimal basis function (OBF) pseudopotential, plane-wave code 10 - for which 
NBSE was originally designed - and in the other case the same quantities are calculated 
from ABINIT. In both calculations twenty valence and conduction bands were included 
with an eight Hartree cutoff criteria for the vectors G. The close quantitative agreement 
between OBF and ABINIT thus serve as a quality check on the various theoretical and 
algorithmic approximations used in our interface. 

Our AI2NBSE interface can also be applied to finite momentum transfer calculations. 
For example, Fig. 2 illustrates the spectrum for a momentum transfer of magnitude q = 0.8 
a^ 1 , where a$ is the Bohr radius. Results from a recent TDDFT calculation are also plotted.— 



IV. OPTICAL SPECTRA OF a-QUARTZ 



The excitonic character of the quartz optical spectra had been recognized even before first- 
principles investigations became possible.— 1 ^ The large band gap and small static dielectric 
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FIG. 2: Calculated €2(q, u>) for Si with impulse magnitude of q = 0.8 clq along the [1, 1, 1] direction. 
AI2NBSE's result (solid line), and for comparison the result from a recent TDDFT calculation 
(dashed line) are plotted. 29 

constant of a-quartz-and the prominent low-energy features in the spectrum-suggested a 
strong excitonic role in the optical properties, which was subsequently confirmed by first- 
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reproduce the measured 



principles calculations.— Our calculations, like those of Ref. 
spectra in detail, with excellent agreement in oscillator strength and calculated feature 
positions corresponding with those measured to within tenths of an eV. The results can be 
seen in Fig. 3. 

The calculations of Ref. |24| were performed for photon polarization within the hexagonal 
plane (i.e., the "ordinary ray"), noting that this is the most commonly measured spectrum. 
This material is however uniaxial, possessing an independent optic axis normal to the hexag- 
onal plane. We have also carried out calculations for photon polarization perpendicular to 
the hexagonal plane (i.e., the "extraordinary ray"). As illustrated in Fig. 3, quartz exhibits 
some anisotropy, but it is less pronounced than for rutile. This result is consistent with 
the known role of the Si04 tetrahedron versus the anisotropic crystal structure, and the 
structural phase insensitivity of the optical functions in quartz. 30 

While birefringence in quartz is a well-studied effect, it is small relative to Ti02 crystals 
discussed in the the next section. Below the interband transition energies, but above the 
lattice response, we find static indices of refraction of 1.52 for the ordinary polarization, and 
1.53 for the extraordinary at 0.7 eV, in precise agreement with the measured values.— 
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FIG. 3: Imaginary part of the dielectric function for SiO"2 with in-plane polarization (top), and out- 
of-plane polarization (bottom). Theoretical spectra from ABINIT and OBF (solid line and dashed 
line) spectra are plotted for both polarizations. The calculated spectra use a 28 Hartree plane-wave 
cutoff and included a 216-point Brillouin zone sampling and the 26 lowest-lying, degenerate bands 
above the Fermi level. The experimental spectrum, reproduced from Ref. 
for the ordinary ray. 



24 is also shown (crosses) 



V. OPTICAL SPECTRA OF RUTILE AND ANATASE 



As one of the simplest transition-metal oxides, TiC>2 exhibits a variety of natural crystal 
structures and presents a fundamentally interesting system for first-principles electronic 
structure methods. This material is also an important component in various ultraviolet 
(UV) applications. It has been demonstrated that ab initio methods can describe various 
physical properties of TiC>2.— ^ There have been a few ab initio studies of UV dielectric 
spectra for rutile and anatase phases. Glassford and Chelikowsky^ reported calculations for 
the rutile phase using a plane-wave pseudopotential approach; Mo and Ching^* used a linear- 
combination-of-atomic orbital method for both phases (and brookite); and Asahi et al. 35 

9 



studied the anatase phase with a linearized augmented plane-wave method. However, none 
of these studies included excitonic effects. A theoretical treatment that includes excitonic 
states may contribute to understanding its optical properties. Toward this end, we present 
calculations of the spectra of Ti02, with polarization dependence, for both the rutile and 
anatase phases using our AI2NBSE interface. 

Each of the two Ti02 phases is tetragonal, and Mo and Ching^ and Fahmi et alM 
have reviewed their structural properties including crystal structures, space groups, and 
differences in bond lengths and angles. These two crystal structures can be considered as 
arrangements of slightly distorted oxygen octahedral elements with a titanium atom at the 
center of each, so that each titanium has an oxygen coordination of six, and each oxygen 
has a titanium coordination of three. The relationship between the two structures has been 
described in terms of varied orientation among the octahedral chains. The two polymorphs 
studied here can be generated with six-atom unit cells corresponding to two TiC>2 units. 
Each unit possesses two inequivalent bonds of "apical" and "equatorial" character, such 
that each titanium sees two apical and four equatorial bonds, while and each oxygen sees 
one apical and two equatorial bonds. 

Structural similarities in rutile and anatase lead to similarities in their electronic struc- 
ture. In Fig. 4 we plot the ABINIT-calculated densities of states for for the two materials. 
The local-density approximation (LDA) gives the uppermost valence-band width of about 
5 eV for anatase and 6 eV for rutile, and this band is regarded to be dominated by O 
2p character. The anatase LDA band gap is greater than that of rutile by about 0.2 eV, 
corresponding with the measured difference.— Because the oxygen bonding environment is 
planar in both phases, there is a decomposition of the O 2p into p^ and p a states.— 1 ^ 

The first conduction band is dominated by Ti 3d character, with the lower half regarded 
as t2 9 -like, and the upper half regarded as e 9 -like.— ^ As can be seen in Fig. 4 and from 
the band structures,— these sub-bands are reasonably resolved energetically. Our calcu- 
lations were performed with semicore states treated as core using Troullier-Martins-type 
pseudopotentials,— ^ and as valence using Teter-type pseudopotentials.— >^ While the semi- 
core states have spatial localizations comparable to oxygen 2p states, their energies are many 
eV below the Fermi level. Our calculations suggest that the optical and UV spectra are not 
highly sensitive to the treatment of Ti semicore states, and that the primary source of the 
discrepancy between theoretical and experimental spectra is the neglect of the excitonic 

10 
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FIG. 4: Density of states p(E) with respect to the Fermi energy Ep for the uppermost valence and 
lowest conduction bands of TiO"2 in the rutile (solid line) and anatase (dashed line) phases. No 
band gap corrections are included. The LDA gap is calculated to be 2.0 eV for rutile and 2.2 eV 
for anatase. 

effects in the low-energy region. 

A. Rutile Phase 

The spectra calculated for rutile with the AI2NBSE interface (Fig. 5) demonstrate sig- 
nificant excitonic effects, as can be seen by comparison with calculations (Fig. 6) which 
neglect electron-hole interactions and are calculated within the RPA. One signature of the 
excitonic effects is the sharp onset in our BSE calculations, which better reproduces exper- 
imental spectra. Our results for the polarization dependence in Ti02 show a much more 
prominent anisotropic optical response than in Si02- This results in strong birefringence, as 
seen in Fig. 7. Our calculated static indices of refraction differ by 0.30, in agreement with 
the observed birefringence. - 

Two pronounced anisotropic features in our calculations are consistent with experimental 
spectra.- 1 ^ First, the primary onset absorption feature assumes a doublet structure for the 
ordinary ray which is not evident for the extraordinary ray. These low-energy excitonic 
features at 4.0 and 5.3 eV are clearly evident in our calculation and are expected to involve 
t2 9 -like final states.— Second, the onset of absorption for the extraordinary ray is significantly 
stronger than for the ordinary. For the extraordinary ray, a singe feature is measured at 4.1 
eV, which is stronger, broader, and more symmetric than for the ordinary ray. However, a 
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FIG. 5: Imaginary part of the dielectric function for rutile with in-plane polarization (top) and 
out-of-plane polarization (bottom). Theoretical from ABINIT and OBF (solid line and dashed 
line) spectra are plotted for both polarizations. The calculated spectra use a 32 Hartree plane-wave 
cutoff and included a 216-point Brillouin zone sampling and the 26 lowest-lying, degenerate bands 
above the Fermi level. For comparison the experimental spectrum from Ref. m (crosses) is also 
plotted. 

third anisotropic feature of our calculations, namely the much stronger absorption peak for 
the extraordinary ray above 6 eV, is not clearly seen in the experimental reflectivity spectra. 
It should be noted that the peak at 7.4 eV is more prominent in the x-ray absorption spectra 
than in the reflectivity data,— ^ and is expected to involve e g -like conduction states.— To 
summarize, our BSE calculation gives a stronger, less structured t 2g -\ike absorption band 
for the extraordinary ray than for the ordinary, in accord with experiment. 



12 




4 6 

Energy (eV) 




4 6 

Energy (eV) 

FIG. 6: Imaginary part of the dielectric function for rutile with in-plane polarization (top) and out- 
of-plane polarization (bottom) calculated without core-hole interactions (Vd = 0). Two spectra 
are plotted corresponding to a treatment of semicore states with Teter-type Ti pseudopotentials 
(solid line), and with Troullier-Martins pseudopotentials where the semicore states are pseudized 
(dashed line). 

B. Anatase Phase 



Theoretically, both the valence-band and conduction-band densities of states are similar 
to that of rutile, as shown in Fig. 4. However, their experimental absorption spectra show 
more variation (Fig. 8.). Strikingly, the measured spectrum of anatase for the ordinary 
ray (in-plane polarized) shows a clear low-energy shoulder below the primary absorption 
feature.— >^ This feature is not seen in rutile, and calculations neglecting the electron-hole 
interactions in anatase do not clearly resolve it.— >^ Aside from this low-energy feature, 
the spectral structure and anisotropy of anatase, both calculated and measured, have some 
similarities to that of rutile. Thus features in the spectra have been analyzed,— 1 ^ and 
conclusions regarding the relevant bands are reminiscent of results for rutile.— For the e 9 -like 
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FIG. 7: Ordinary (dashed line) and extraordinary (solid line) indices of refraction n of rutile-TiO"2 
as a function of photon energy. At 589 nm, which is below the band gap but well above the lattice 
resonances, the calculated difference between the two indices is 0.30 while the measured value is 
0.29.-2 

absorption band, the experimental spectrum for anatase does show a stronger absorption for 
the extraordinary ray, 43 similar to what is calculated in both the rutile and anatase phases, 
but not seen in experimental rutile spectrum. 

C. Interpretation of Polarization Dependence 

Molecular orbital interpretations of the Ti02 spectra have been reported widely 4 ^ 4 - and 
provide a qualitative picture of the electronic structure in this system. Such analyses illus- 
trate the hybridization between the Ti 3d and O 2p states. However, immediately above and 
below the Fermi level, the hybridization has been shown to be weak.— ^ The calculations of 
Sorantin and Schwartz 35 demonstrate that the valence band for rutile just below the Fermi 
level is primarily 7r-bonding O 2p character, while lower in energy is primarily cr-bonding O 
2p character. Just above the Fermi level, the states are, in energetic order, primarily Ti £ 2 g 
and e g . 

A molecular-orbital analysis, considered in conjunction with measured x-ray spectra, al- 
lowed an empirical identification of features in the optical spectra of Harbeke and Cardona 9 
with specific electronic transitions.— Although selection rules, inferred from the approximate 
atomic states, were invoked in that work, no effort was made to analyze the polarization de- 
pendence of the selection rules. Some of the most prominent features could not be identified 



14 




4 6 

Energy (eV) 




4 6 

Energy (eV) 



FIG. 8: Imaginary part of the dielectric function for anatase with in-plane polarization (top) and 
out-of-plane polarization (bottom). Theoretical spectra from ABINIT and OBF (solid and dashed 
lines) are plotted for both polarizations. The calculated spectra use a 32 Hartree plane- wave cutoff 
and include a 216-point Brillouin zone sampling and the 26 lowest-lying, degenerate bands above 
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the Fermi level. For comparison the experimental spectrum reproduced from Ref. 9j (crosses) also 
plotted. 

at all, and were attributed to the unaccounted excitonic effects. Assignments of the same 
experimental features were made, for ordinary and extraordinary rays independently, from 
a band-structure perspective by Glassford and Chelikowsky. 33 But the calculated spectra 
again did not treat the excitonic interactions, and hence did not reproduce the experimental 
features to the level of agreement reported here. Also, local electronic structure was not 
emphasized in the peak assignment analysis. 

Perhaps the most conspicuous anisotropic characteristic in the rutile spectra is the 
stronger onset absorption for the extraordinary ray. We observe that the first eV below 
the Fermi level is dominated by the O 2p n states, and further, that the oxygen bonding 
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planes are all defined by in- hexagonal-plane normals. This implies that, for some in-plane 
polarization, there is a single dipole allowed channel from one of the two 7r-state orienta- 
tions (corresponding to the two oxygen-bonding planes), and the three ti g states, while for 
the extraordinary polarization, there are two dipole-allowed channels. This is consistent 
with the calculated and measured stronger threshold resonance for the extraordinary ray. 
In connection with this point, we emphasize that while the t 2g and e g symmetry labels are 
approximate because of angular and bond-length distortions within the octahedra, the se- 
lection rules above are unmodified from the above even when the local, D 2 h symmetry group 
is considered, while the actual crystallographic group in rutile is of even higher symmetry. 

VI. SUMMARY AND CONCLUSIONS 

In summary, we have developed a first principles approach for calculations of various op- 
tical spectra, including finite momentum transfer in crystals. The method combines ground 
state electronic structure calculations from ABINIT with BSE calculations from the NBSE 
code. The method is tested on bulk Si, yielding results in good agreement with other meth- 
ods. Calculations are presented for the macroscopic dielectric spectra and its orientation 
dependence in a-quartz, rutile Ti02, and anatase TiC^. Our quartz spectra for the ordinary 
ray reproduce the strong excitonic character and are in good agreement with experiment 
and earlier theoretical work. The anisotropy of the rutile and anatase phases of TiC>2 are 
more significant than for quartz. The static indices of refraction for the two polarizations 
of rutile differ by more than 10%, in agreement with experiment. Also our calculated ab- 
sorption at low energies reproduces experiment better than previous theoretical results that 
neglect excitonic effects. In particular our calculations reveal an additional low-energy fea- 
ture in anatase also found in experiment. Also we are able to interpret the anisotropy in 
the threshold behavior for rutile in terms of n to £29 selection rules. 
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APPENDIX: AI2NBSE INTERFACE 

In order to calculate the optical spectra, the AI2NBSE interface first obtains Kofm-Sham 
energies and wave functions from the self-consistent ground-state electronic structure code 
ABINIT. The single-particle eigenenergies are modified to include self-energy corrections 
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according to the GW approximation. Subsequently the interface constructs several quantities 
needed for the NBSE calculations. These include current-operator matrix elements between 
Kohn-Sham states and the ground-state charge density for calculating Hybertsen-Levine- 
Louie screening. Thus our interface does not take advantage of ABINIT's dielectric function 
capability. 

Typical AI2NBSE calculations are divided into four modular stages and require a single 
input file which contains all parameters needed to define both the system and its ground- 
state and excited state one-electron properties. Briefly the modular operations are as follows 
after the input file is read and stored: 

1) ABINIT calculation: ABINIT inputs are generated and ABINIT is run. These cal- 
culations supply the Kohn-Sham eigenvalues e„k and eigenfunctions -0„k both for the oc- 
cupied and unoccupied states. Currently the interface uses a standard serial distribution 
of ABINIT. However, calculations for large systems can still be performed by means of 
consecutive runs which are automatically set up by the interface, with minimal input or 
intervention from a user. 

2) Density components: Fourier components of the ground-state density p(G) are generated 
using the eigenfunctions from stage 1. 

3) Dipole matrix elements: Dipole matrix elements in Eq. (8) are calculated, and the eigen- 
functions are converted to the format used by NBSE. 

4) NBSE calculation: All quantities required for the Bethe-Salpeter calculation are collected 
and NBSE is run. The output includes various optical spectra and optical constants. 
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